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Abstract 

Boundary conditions may change the phase diagram of non-equihbrium sta- 
tistical systems Uke the one-dimensional asymmetric simple exclusion process 
with and without particle number conservation. Using the quantum Hamilto- 
nian approach, the model is mapped onto an XXZ quantum chain and solved 
using the Bethe ansatz. This system is related to a two-dimensional vertex 
model in thermal equilibrium. The phase transition caused by a point-like 
boundary defect in the dynamics of the one-dimensional exclusion model is 
in the same universality class as a continous (bulk) phase transition of the 
two-dimensional vertex model caused by a line defect at its boundary. 
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It is generally held that for statistical systems the influence of boundary conditions 
should be negligible, at least when one is interested in the bulk system only. This 
belief is indeed supported for equilibrium systems by abundant evidence coming from 
analytical or numerical studies. As we are going to show in this letter, however, this 
picture is far from being generally valid. Boundary conditions may indeed be crucial 
for the bulk properties and may even cause continous phase transitions. 

It has already been realized that a change in boundary conditions (equivalent to 
some localized point defect in a system on a ring) may cause various kinds of phase 
transitions in the static properties of one- dimensional non- equilibrium systems [|l], 0. 
Here we show that boundary terms may also induce phase transitions in the dynamics 
of these systems. Such phase transitions are then shown to correspond to bulk phase 
transitions of two-dimensional equilibrium systems caused by boundary terms (i.e., 
line defects). 

The paradigmatic example we are going to study is the asymmetric simple exclu- 
sion process, see . In its simplest version, without particle creation or annihilation, 
it is described by particles of a single species A moving on a lattice. A given site j can 
be occupied or empty at an instant of time t. A particle at site j for time t may hop 
at time t + 1 to its right neighbor with rate (1 + e)/2 and to its left neighbor with rate 
(1 — e)/2, if the final site is empty. This simple model appears in a large variety of 
contexts. It has been argued to be in the same universality class as the noisy Burger's 
equation 0]. This in turn is a one-dimensional version of the incompressible Navier- 
Stokes equation or else can be regarded as the one-dimensional Kardar-Zhang-Parisi 
equation describing the shape fluctuation in various growth models. 

Using a master equation approach, the probability distribution function P{{/3}]t) 
is obtained by solving dtP = —HP, where /5(t) is a configuration of occupied and 
empty sites and the quantum Hamiltonian reads 



where L is the number of sites and a^'^'^ are Pauli matrices. The ground state energy 
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(1) 



vanishes as follows from probability conservation in the master equation. Obviously H 
commutes with the particle number operator = ^2^=1 where we have introduced 
the particle projector rij = |(1 — cr|). The particle density is then p = N/L. Using 
periodic boundary conditions = crl'^'^, Gwa and Spohn [Q then show using 

Bethe ansatz techniques for e = 1 and p = 1/2 that the (real part of the) energies 
El for L large scale as 

El ~ L-'/' (2) 

On the other hand, the Hamiltonian had also been derived considering asymmetric 
lattice diffusion for the case of free boundary conditions where the particles are not 
allowed to move beyond the boundary (impenetrable walls). Then it reads 

L-l 



H' = y 



(3) 



where 



This is the well-known XXZ quantum chain which is symmetric under the quantum 
group UqSU{2) as shown by Pasquier and Saleur We stress that the quantum 
Hamiltonian of the diffusion process on a periodic lattice cannot be obtained by 
simply taking periodic boundary conditions in eq. (0). Using the UqSU{2) symmetry 
it is easy to show that for e 7^ the energies for L large become 

~ 1 - (5) 

in each sector with particles, see [|, 0, §]. 

The apparent inconsistency of eqs. (0,^) will be explained below. Our analysis 
applies to any value of e, generalising the approach of [Q. We note that since the 
Hamiltonian of eq. (|l]) is gapless, the approach of time-dependent mean values towards 
their equilibrium value may be according to a power law, while with the Hamiltonian 
of eq. (§) the relaxation will be always exponential. 
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To understand this observation, we begin by rewriting eq. (|I]). Introducing raising 
and lowering operators cr^ 



|((T^' ± icpy), we have 



H 



1 +e 



E 



+ - 



a, a.+i + 



2fl 



— = — y 



+ - ,-!- + , Q + Q 



-1 



^2 



1 



(6) 



where the second of the eqs. (^) has been used. Next, consider the non-singular 
operator 



f/ = exp ngJ2j(^j 



±TT-1 



UafU 



(7) 



If we now choose q = e^'^^, we obtain 



H" = UHU 



-1 



2(g + g- 

which is indeed (almost) the Hamiltonian H' . The distinction comes from the bound- 
ary conditions. The surface field (g — q^^){al — cr£)/8A is absent in H" and one has 



a 



L+i = q^ <y 
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cr 



L+l 



(9) 



which make H" non-hermitian, as H already is. It is these unusal boundary conditions 
which give rise to the different properties of H and if', as we now show. 

First we note that the ground state of the system with free boundary conditions 
as well as that of the model (||) with twisted boundary conditions (|) is (L + l)-times 
degenerate and has energy 0, independent of L. In what follows, we assume L to 
be even. In each sector with fixed particle number N the lowest energy is 0. This 
can be shown by either solving the Bethe ansatz equations or by using the symmetry 
relations with the generators of the quantum algebra |§, |^. 

The calculation of the spectrum proceeds via the Bethe ansatz. Indeed, the XXZ 



chain with non-periodic boundary conditions was already considered in |T0| and we 
merely have to adapt their results to the problem at hand. To begin with, we consider 
the sector with = 1 particle. Then the energies are in this sector 



E = 1- A-'cosO , e = 2-n[ig + - 



n 



(10) 



where n is an integer from the set {0, ±1, . . . , ±(L/2 — 1), L/2}. For L large, the 
energies become 

cos(27r {ig + l)) 

h/ = i 

cosh(27r5f) 



~ 27r^ (^-j +27r?tanh(27r^)- + ... (11) 

and we note that the mass term indeed cancels. We also observe that the real part 
"RE ~ L^^ and the imaginary part '^E ~ L^^. We return to an explanation of this 
below. 



Next, we take the sector N = 2. From the Bethe ansatz |TD[ we have 



E={2A- cosO - cos d')/\-^ (12) 
eL-2TTigL = 2nI -Q{e,e') , O'L - 2nigL = 2nr - Q{e',e) 

where /, /' are distinct half-integers from the set {±|, ±|, . . . , and 

9 — 9')^ 

e{9, 9') = 2 arctan I . / ^ I (13) 



A sin 


{1(9-9')) 




cos 


(1(^ + ^0) 


— A cos 


{W-9')) 



Define 9 = 9 — 2iTig, 9' = 9' — 2-Kig. Then, for small values of the arguments 

z coth(27r^) ■ I — (14) 

which is of order unity. Therefore 

9 = ^a , 9' = ^a' (15) 

where a, a' are of order 0{1). It follows that the same cancellation as observed in the 
sector = 1 also takes place here and also that the observed scaling of the energies 
does not change. Finally, for arbitrary 

N 



E = A-^ (^NA - J2^cos 9r}j 



(16) 



with 

N 

9mL - 2mgL = 27r/^ - J] Q{9^, 9^) (17) 

ra=l 



and the same argument can be repeated. It follows that for any value of e, the 
spectrum is massless. There is a distinction, however, between finite particle densities 
p = N/L = 0{1) and small densities p = 0{L^^). As shown above, in the latter case 
the real part of the energy gaps scales as L^"^, while for p = 1/2 it vanishes like 
(when e = 1) [^]. The way how this result was achieved suggests that this scaling 
behaviour is characteristic for all finite densities. 

After this analysis of the spectra of H' and H" (or H, as the spectra of H and H" 
are identical) we turn to an interpretation of our results. The iV-particle ground state 
|A^) is stationary with respect to the stochastic process defined by the Hamiltonian 
and we shall refer to it as the steady state of the system. Average values (C) of 
operators C in some A^-particle probability distribution \Pn) are defined as (C) = 
{N\C\P]\f). Of particular interest are the n-point correlation functions in the steady 
state 

G{xn, . . . ; xi, h- 0, 0) = {N\n,„e-''^'--'-'^ . . . n^.e-'^'^uolN) (18) 

and their connected counterparts Gc defined by replacing the operators by {rix — p) 
in dg). 

We discuss the system dynamics by studying the dynamic structure function 
SN{k,t) which is the Fourier transform of the connected two-point correlation func- 
tion in the A^-particle sector. In the case of asymmetric diffusion with free boundary 
conditions (|^) the correlation function decays exponentially in time with finite relax- 
ation time r = A/ (A — 1), see eq. (||). 

SN{k,t) oc f{k)e-'/^ (19) 

The relaxation time does not depend on the density p. In the first model, however, 
corresponding to asymmetric diffusion with periodic boundary conditions, the relax- 
ation time diverges. For small densities (finite number of particles) we read from 
eq. (|1T]) r ~ while the result obtained by Gwa and Spohn translates into r ~ 

^3/2 

for finite densities. For the dynamic structure function this suggests the scaling form 
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S{k, t) oc e-''^^-^P^'''h{0H) (20) 

for finite densities while computing the exact structure function for one particle from 
eq. (pup suggests 

S{k,t) oce-''%-'^'''' (21) 

for any finite number of particles. The phase factors exp (— ^e(l — 2p)t) arise from 
the imaginary part of the energy gaps and refiect the steady state current of particles 
moving around the ring (recall tanh{2TTg) = — e). The scaling functions h{k^^'^t) and 
exp(— /c^t/2) have their origin in the diffusive nature of the process. From these 
expressions one can read the dynamic exponent of the process: it is z = 2 for small 
densities and z = 3/2 for finite densities. It is interesting to note the consistency 
with the results for the diffusion constant D of a single tagged particle in the fully 
asymmetric exclusion model studied by Derrida et al. [Tl|]. These authors consider an 
exclusion model of the kind discussed here containing one particle which has the same 
dynamics as all the other (indistinguishable) particles, but is tagged. (Effectively, this 
a model with two types of particles, A and B, containing A^i particles of type A and 
Nb = 1 particle of type B all moving with the same probability to the right if the 
neighbouring site was empty.) They follow the motion of the tagged particle and 
compute its diffusion constant D. They find that D is finite for a finite number 
of untagged particles, corresponding to an dynamical exponent z = 2. On the other 
hand D diverges proportional to L^/^ for finite densities which corresponds to z = 3/2. 

The discussion of behaviour of the dynamic structure function for free boundary 
conditions (|^) on the one hand and for twisted (periodic) boundary conditions ( pO|) 



and (21) on the other hand shows that a phase transition takes place when changing 
the boundary conditions of the model and elucidates its effect on the dynamics of 
the system. Now we show that this corresponds to a bulk phase transition of a two- 
dimensional six- vertex model. The Hamiltonians if^ and H'l of the model with L sites 
can be derived from the transfer matrix Tl of the six-vertex model with a defect line 
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of the type given in eq. (j^) located between sites L and 1 (which we call the boundary) 
0, ^ . Its partition function Z in thermal equilibrium on a L x M lattice is given by 
the trace of the M*^ power of Tl. In the thermodynamic limit L, M cxd one obtains 
Z = lim^ M->oo Tr [exp (— if^M)]. From the vanishing of the energy gaps discussed 
above (which corresponds to the appearance of an infinite degeneracy of the ground 
state of H) we conclude that changing the boundary conditions causes a continous 
phase transition in the two-dimensional six-vertex model in thermal equilibrium. 

So far we have studied asymmetric diffusion with a conserved number of particles. 
We briefly show that a boundary induced phase transition occurs also in an asym- 
metric diffusion model with pair annihilation A + A —>■ We study the Hamiltonian 



'L-l 



2(g + g-i 

with 



H = -7;7Z^—rAT.h, + bh,] (22) 



/i . = ga+(T^.+i + q V^. a+^, + " 2) + Q^j + Q^'^U ■ (23) 

The boundary conditions are defined by the parameter b. If b = one has free 
boundary conditions (we call this Hamiltonian when a distinction w.r.t. the 
boundary conditions is necessary) while 6 = 1 correspond to periodic boundary con- 
ditions (denoted H^). H describes a process where particles in a pair of sites hop 
with probability (1 + e)/2 to the right and probability (1 — e)/2 to the left if the 
respective sites are empty. Pairs of particles are always annihilated. H has a twofold 
degenerate steady state with energy 0. These are the ferromagnetic ground state 
containing no particles (all spins up) and the one-particle state (one spin down) in 
which each possible position of the particle is equally probable. 

The operator U eq. (^ transforms into the free fermion Hamiltonian studied 
in which describes a transition of the Pokrovsky-Talapov type [Q. The energy 
gaps are all finite and of the form 

~ iV f 1 ^ 1 (24) 



q + q 



for large L, which imphes an exponentially slow approach to the stationary state at 
late times. 

For periodic boundary conditions is transformed by U into 



1 



H'" = - ^ ^ , V 

2(g + g-^),tl 



(25) 

with the twisted boundary conditions eq. (P). Following the argument of the 
spectrum of H^' is seen to be independent of the annihilation term a^o^j^^ and it is 
therefore equal to the spectrum of 



1 ^ 



(26) 

which commutes with the particle number operator and describes non-interacting 
fermions on a ring with twisted boundary conditions eq. (^). The eigenvalues of H^' 
are easily found by a Jordan- Wigner and Fourier transformation and are of the form 

E"" = y(l ^ cos d\ (27) 

similar to eq. (|T6|) but with 

= 271 [ig + mjL) (28) 

where the are pairwise distinct integers < mj < L — 1 for odd and half- 
integers I < ?77.j < L — I for even. The low lying energy gaps (A^ finite) vanish in 
the thermodynamic limit as '^E^ ~ L^^ and '^E^ ~ L~^. Here, as opposed to the 
model with free boundary conditions, finite-density states decay with a relaxation 
time of order 1, while for low-density states (finite A^) the relaxation time diverges 
proportional to L^. 

The Hamiltonian ( ^2]) is related to a 7- vertex model with a boundary defect, 
see As in the 6- vertex model discussed above, changing the boundary condition 
induces a bulk phase transition. The same mechanism should also apply to more 
general reaction-diffusion problems, for example two-particle processes with a reaction 



y4 + i? — > ||T3l, or even extensions of the Hubbard model |jT4|, where the same 



quantum chains as considered here reappear 
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To summarize, we have shown that boundary conditions can have a major in- 
fluence on the phase diagram of certain statistical models. Since the models con- 
sidered here are merely prototypes of much more general ones, we expect that the 
phenomenon found is generic. 
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